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1. Introduction 

Despite progress since the early 1980's, when it was realised that perturbative tech- 
niques fail to capture its critical behavior, the Random Field Ising Model (RFIM)[|l| is 
still in need of further illumination. The cause of the difficulty in dealing with this type of 
disorder is that the energy landscape is complicated, with many local minima that pertu- 
bation theory fails to take account of. Mean field theory provides a simple insight into the 
difficulty: at sufficiently low temperature the mean field equations have many solutions. 
The intuitively obvious way of defining mean field theory would be to weigh these solutions 
according to their Boltzmann factors. This weighted mean field theory, or to be precise, 
the theory with the weights e~^^ which include an entropy factor, is the subject of this 
paper. In this context the failure of the pertubative approach comes about because, as 
is clear in the supersymmetric formulation^, the prescription for the weights does not 
contain a Boltzmann factor. 

In this report we solve the mean field equations numerically and directly construct 
the weighted mean field theory. Our main interest is to see how critical behavior can 
arise at the ferromagnetic transition. Guagnelli et aZ. have analysed a truncation of 
the theory by only considering the maximal solutions: those with maximum or minimum 
magnetisation. They found that critical behavior does not arise in that case. Here we shall 
enlarge on their work by considering a much larger number of solutions and shall see how 
divergences in thermodynamic quantities can arise. The number of solutions included in 
the sum is clearly important in this approach and we shall discuss the issue in detail. 

Replica theory is the natural setting for dealing with the difficulty mentioned above 
and although not yet successful in treating the critical behavior, it is helpful to keep 
the resulting scenario in mind. The replica symmetric solution becomes unstable as 
one reduces the temperature below Trsb, at which point the correlation length remains 
finite0. The critical temperature Tc at which ferromagnetic order arises is at a lower 
temperature Tc < TfisB- Finally one expects replica symmetry to be restored at an even 
lower temperature. This scenario has clear analogues in our mean field treatment. 

We start in sections 2 and 3 by discussing the consequences of weighted mean field 
theory and describing the technique of numerically searching for solutions. Section 4 is 
concerned with the general properties and number of solutions we find. Analysis of the 
correlation functions and the emergence of critical behavior is the subject of section 5. 
Finally we present a short conclusion discussing the numerical values of the exponents we 
determine. 
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2. The RFIM and its Weighted Mean Field Theory 

The system consists of Ising spins subject to quenched random fields hf, 



H = -Y1 S,Sj-Y,h^S,. (2.1) 

{i,3> 

The effect of the fields is to destroy the tendency to long range order, and the original 
argument of Imry and Ma0 based on the energy balance for domain formation gives the 
lower critical dimension as died = 2. Perturbative methods, which lead to the phenomenon 
of dimensional reduction [^|^, would instead predict died = 3. However in three dimensions 
it has been shown rigorously that long range order prevails at low enough temperature]^. 
This case is most interesting and we shall restrict ourselves to studying the system in three 
dimensions. 

In the work described here the fields are taken to be ±|/i| with equal probability. Ac- 
cording to arguments of Aharony ||] , and consistent with the zero temperature simulations 



of Ogielski [|T0| , a bimodal probability distribution of this sort causes the phase transition 
to be first order when \h\ is very large. As in we work with \h\ = 1.5, which is small 
enough to give a continuous transition, yet large enough to shift Tc substantially from its 
value without disorder, thus revealing a range of non-trivial critical behavior before the 
crossover to pure behavior at higher temperature. 

For a given realisation of the random fields, the mean field equations for the system 

are: 

rui = tanh {(3{Dmi + hi)) , (2.2) 

where is the local magnetisation and Drrii indicates a sum of the magnetisations over 
the nearest neighbours to site i. As has already been indicated, at low temperature there 
can be many solutions to this equation. We shall denote the solutions by mf, and in general 
shall use the superscript a to denote a quantity, such as the free energy F", calculated for 
that solution. 

= V (^F" - , (2.3) 

where the energy F", and entropy S", of the solution are given by the expressions, 

(2.4) 
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A convenient way of organising the weighted mean field theory is in terms of a mean 
field partition function defined as: 



J2e-^^\ (2.5) 



This way of writing the theory simply encodes the intuitive definitions of quantities, such 
as the energy, as Boltzmann weighted sums over the solutions. 

^ V dp V ^ dp"-^ ' Zmf ^ Zmf ^ 

a a a 

where we have defined the weights: = e~^^ /Zmf- Note that dF°'/d/3 is simple only 
because of the constraints dF°'/dmf = 0, which are nothing other than the mean field 
equations 

The average over the random fields is performed as the last step and is denoted by an 



over bar, for example {E). 

Critical behavior becomes apparent through study of the correlation functions. In 
the RFIM the correlators (SiSj) and {Si){Sj) are both more singular in momentum space 
than the connected correlator {SiSj)c = {SiSj) — {Si){Sj). Although different from the 
field theory usage, it seems common in random field systems to call the correlator {Si){Sj) 
"disconnected" and we shall follow this convention. In section 5 we shall determine the 
exponents fj and 77 associated with the disconnected and connected correlators. 

Contrary to ordinary mean field theory, the probability distribution implied by Zmf 
is not factorised. The fiuctuation dissipation theorem (FDT) is therefore not needed to 
calculate the correlation functions. Using Zmf we find: 

{Sr){S,) = {J2^-mT){J2^^rn]) (2.7) 

a 7 

{SrSj)c = (5^«;"m>-) - iJ2w"mf){Y,^'ym]) + ^ J^^^^', (2-8) 

a a 7 a 

where gfj is defined as dmf /dhj and is the usual term arising from the FDT. 
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3. Solving the Mean Field Equations 

3.1. Iteration Technique 

Starting from some seed configuration rrif \ the mean field equations ( p.2|) are solved 



by iteration 



(t+i) 



tanh (^/3(Dmf ^ + /i,)) • (3.1) 



This is implemented as an efficient code running on the APE parallel processor. We 
insist on strict convergence requirements, that for each solution ^^^oanA^i^^^ —rPbp)'^ < 
2.10"^'^. Where the m^'s are represented to ffoat accuracy and there are 128 APE boards to 
cover the complete 32'^ lattice. With this requirement we find no difficulty in distinguishing 
solutions. The criterion for saying that two solutions are the same, that the maximal site 
difference | mf — m] \ , be less than some cutoff, leads to the same identification of solutions 
for a wide range of cutoff values. We have finally chosen this cutoff to be 10"*^. Although 
the majority of seeds quickly converge, for some temperatures and some realisations of the 
random field we found that the convergence time was unacceptably large. A maximum 
number of 15000 iteration steps has been imposed, leading to rejection of about 1.4% of the 
potential solutions. We have observed a phenomenon known as "funneling" in which the 
configurations arising from different seeds rapidly converge to very similar configurations 



which then follow the same path though configuration space before finally converging |T2 
The iteration technique may lead to solutions that are maxima of the free energy besides 
minima. Although we do not see such solutions in the cases described below, we have no 
simple way of rejecting them and they are implicitly included in the weighted sum. 

3.2. Maximal Solutions 



The iteration (|3.1|) has the property that if one can assign an ordering to two con- 



figurations nia and m^: '>TT'a\ ^ ^1*1 sites i; then this ordering is preserved: 

> Consequently one can identify two special, maximal solutions, m_|_ 

and m_, that arise from the seeds, m^^j = 1 and m^^j = — 1, and that bound any other 
solution since m^] > m^^y seed i — "^-1- These were the solutions analysed in [Q. 

The existence of maximal solutions provides an accurate means of searching for the 
temperature (which by analogy with the replica analysis, we denote Trsb) at which the 
mean field equations start to have more than one solution. Below this temperature m+ and 
m_ differ. We have found Tusb for series of random fields both on 32^ and 64'^ lattices. 
At the larger size the peak of the distribution of Trsb's moves to higher temperature 
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while the width shrinks - supporting the results from replica theory that T^sB is well 
defined and separated from Tq in the thermodynamic limit. In about 10 percent of the 
cases we examined we observed two or more values of Tjisb- That is: at high temperature 
there is a unique solution, then as the temperature is reduced we first find more than 
one solution, then an interval where again only one solution exists, before finally reaching 
another low temperature region with many solutions. This effect is reminiscent of some of 



the observations of jumping made by Sourlas at zero temperature [13 



3.3. Seed Strategies 

We have considered several strategies for the choice of seeds. The most efficient 
strategy we found was to consider chequerboard seeds: we divide the 32"^ cube up into 8 
subcubes of size 16"^, and colour each one independently with +1 or —1. The leads to 256 
seeds including the maximal ones. Further subdivision was thought to be impractical, so 
further sets of 256 seeds were generated by adding independent random pertubations to 
each site of the chequer seeds. 



4. Number and Properties of Solutions 

The number of solutions found and included in the weighted mean field theory is an 
important parameter. We start this section by displaying the average number and the 
effective number of solutions we find in our work. The main justification however, for 
curtailing the search for further solutions at the point we choose, is that the quantities 
measured from the correlators show very little change on increasing the number of solutions 
beyond those found with chequerboard seeds alone. 

The results we shall present are based on data for the 32'^ periodic lattice for a set of 
150 different magnetic field samples. We have calculated quantities based on the maximal 
solutions alone, and the solutions obtained from chequerboard seeds. Besides these, in the 
temperature region identified as most interesting, we have performed longer runs including 
chequerboard followed by two sets of chequerboard plus random seeds for the first 100 of 
the magnetic field samples. The temperature ranges covered by each of these three data 
sets are most clearly seen by looking ahead to the ranges of the three curves of fig. 2. 

In all cases the temperature is measured in units of the mean field critical temperature 
for the pure model T/Tc pure = T/Q. 
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Quantities such as the Energy, Entropy and Free Energy vary smoothly with tem- 
perature and the curves show httle variation when the number of solutions included is 
increased. The average squared site magnetisation is also a smooth function, but the av- 
erage magnetisation itself suffers large fluctuations. The specific heat starts to develop a 
cusp at the temperatures we shall later identify as the ferromagnetic phase transition. 

4.1. Number of solutions 

The number of distinct solutions found by iteration is shown in fig. |l|. Those solu- 
tions with completely negligible weight (F" — Fmin > 80) have been dropped. The three 
sets of points are from successively increasing numbers of iteration trials. The lowest set 
(squares) is for the 256 chequerboard seeds while the upper two sets (triangles) corre- 
spond respectively to the addition of one and two sets of random pertubations on top of 
chequerboard. 

Already from this figure it is clear that there is a law of diminishing returns relating 
the number of solutions found to the number of seeds iterated. From more detailed studies 
at a fixed low temperature, we see that after a steep initial rise, the number of solutions 
only increases slowly with the amount of effort. We have determined to stop searching for 
solutions at the point determined by the shoulder of this curve, which in practice is after 
the chequerboard and two sets of chequerboard plus random seeds. At higher temperatures 
a less though search is sufficient. 

A feeling for the degree to which extra solutions are important comes from taking 
account of their weights, and a convenient quantity is defined by: 



The inverse gives the effective number of solutions that contribute. A plot of is 
shown in fig. |^ for maximal, chequerboard and chequerboard plus random seeds. At high 
temperatures = 1 since only one solution exists, whereas the rise at low temperatures 
is due to the increasing dominance of certain solutions. This increasing dominance can 
also be seen in the density of solutions in free energy which rises less steeply at lower 
temperatures. The fact that the effective number of solutions grows more slowly than the 
total number of solutions gives us our confidence that we are finding all the important 
solutions. 

A histogram of the overlaps between different solutions displays a peak at small overlap 
with a long tail. When the histogram is normalised by the product of the weights of the 
solutions the peak moves to smaller overlap. 




(4.1) 
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4-. 2. Form of solutions 

There is a regime of temperature where it is possible to get some intuition into the 
form of the solutions. Generally the solutions are complicated: at high temperature rrii 
is small but follows the field rrii ~ (3hi + whereas at low temperatures there is 

an overall magnetisation and sometimes reversed field domains. However, in the region 
just below TjiSB where there are only a small number of solutions, it is possible to see 
that the differences between solutions are local. In fig. |^ a three dimensional picture of 
the difference in magnetisation between the maximal solutions (note that this is always 
positive) is shown by the constant difference surfaces. The dark region has small difference. 
The notable feature is that the regions of appreciable difference do not touch each other and 
are separated by a sea where the solutions are almost identical. This suggests that when 
the difference between maximal solutions has local regions, each can be independently 
switched to either of two configurations, and that there should be a total of 2^ solutions. 
In fact, in the example corresponding to fig. |^, we do find a total of 8 solutions, and this 
method has been used to test the efficiency of different seed strategies. These solutions are 
found where the regions are local and do not infiuence each other. At such temperatures 
we must expect a constant density of such local differences between maximal solutions, and 
thus a total number of solutions that grows exponentially with the volume of the system. 
This picture fails at lower temperatures where the difference between maximal solutions is 
no longer local. 

As was the case for the maximal solutions 0, in any particular solution correlation 
functions are rather rough with finite correlation length. 

4-3. Below the transition 

The reader will have noticed that our definition of mean field theory as a sum over 
all states is at variance with our usual understanding of the the ferromagnetic phase 
transition. One might try to modify the definition at low temperature by only including 
solutions with one sign of the magnetisation, however the region in the vicinity of the 
transition will remain unclear. It is possible to visualise how the transition takes place 
from the lists of solutions along with their magnetisations and weights. In the top row of 
fig. H the number of solutions is shown as a histogram against their magnetisation. The 
width of the distribution increases at lower temperature, but there are always solutions of 
small magnetisation. In the lower line of plots in fig. ^ the summed weights rather than the 
number of solutions are shown. It is clear that at low temperature the small magnetisation 
solutions have low weight and are unimportant, and that the space of solutions is divided 
into two significant groups with opposite magnetisation. 
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5. Critical Behavior 



In order to study critical behavior of the theory we look at the susceptibilities and the 
quantities that arise from analysing the correlation functions. Besides investigating the 
divergence of the correlation length, we expect two independent exponents defined by the 
behavior at the critical point. In three dimensional real space the disconnected case goes 
as ~ r^~^ whereas the connected case dies more quickly as ~ r~^~'^. 

We always consider plane-plane correlators in one spatial direction with the transverse 
momenta set to zero: G{x,ky = 0, /c^ = 0). From the solutions mf two different lattice 
correlation functions, C*-^-* and C^'^\ can be defined using sums over transverse planes. For 
example in the x direction: 

i^"(^) = 7VE^"(^'2/'^)- (5.1) 

f y,z 

For C^^\ which is the plane-plane version of the disconnected correlator (|2.7|) , the weighted 
sum over solutions is, 

(5.2) 

CW(x) = 5^u;"CW"(x). 

a 

Instead of performing the sums in the order given above, we can first weight the plane 
sums before combining them, 

P.(x) = 5^«;"P,"(x) 

1 (5.3) 
Ci^H^) = -^^Px(a;o)J"x(xo + x). 

This procedure yields the plane-plane version C*^^-* of the first term in the connected cor- 
relator {SiSj)c, (|2.8|). An average over the principal directions and over the random fields 
is made before obtaining the final results for and \ The plane-plane connected 
correlator is given by C*-^-* — C*^^-* + C^^^'^\ in which the FDT term cannot be determined 
directly from the solutions mf and a further iteration is required. 
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5.1. Disconnected Correlator, {Si){Sj) 

Scaling arguments suggest that away from the critical point the 3 dimensional discon- 
nected correlator behaves as G{r) = r^~^ f{r /^). The scaling function we shall use to fit 



that is motivated by zeroth order pertubation theory. This choice is recommended by its 
simplicity and it is frequently employed for fitting experimental data. The disadvantage 
of this choice of scaling function /(r/^) = {^/r)^~'^e~^^^ , is that the critical limit, ^ oo, 
amounts to 77 = 1 behavior. Although 77 is in fact close to 1, we have explicitly checked 
that the other natural choice of scaling function f{r/^) = e~'^^^, gives similar results. 
The plane-plane form of our fitting function is: 



In contrast to the lack of sensitivity to the precise form of scaling function, it is 
important to make the periodic modification correctly. Since the correlations fall off so 
slowly, only keeping the leading two terms in the periodic sum leads to significant errors. 
We use the complete form: 



In all cases, even well away from the critical point, the fits are extremely good. Al- 
though we work exclusively above the critical temperature a constant term B has been 
included. This is done in order to take account of finite size effects and is also necessary 
to avoid discontinuities between data sets from different number of seeds. In the infinite 
system we would expect B to be the mean site magnetisation squared m^; we however 
always observe it to be larger, approaching m^ at low temperature. We shall regard B 
becoming different from zero as a signal of finite size effects being important. 

In fig. ^ and fig. ^ the constant B and the correlation length obtained from the fitting 
are plotted against temperature. The jacknife technique is used for the errors, but note 
that they are correlated since the same set of random fields is used at each temperature. 
The plots are displayed for the maximal and for the chequerboard solutions. These sets 
of data fiow smoothly into each other as do the points based on the solutions arising from 



C*-^^ and define the correlation length arises from a Lorentzian squared propagator, a form 




(5.4) 




(5.5) 
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chequerboard followed by chequerboard plus random seeds. In fact these later points lie 
well within the errors of the chequerboard points and are not shown. It is generally the 
case that extra solutions beyond chequerboard do not change any of our analysis of the 
disconnected correlator. For each of the fitting parameters the quantity a;*/3(a;^)^, remains 
close to 1 (its value for a gaussian distribution), except for the measurement of the constant 
term B at high temperatures where it becomes as large as 2. 

Consider fig. |] for the magnetisation squared as measured by the constant term B of 
the correlator. Significant deviations between the points calculated with many solutions 
and those calculated with only maximal solutions become apparent below the temperature 
at which the maximal solution B becomes non-zero (T ~ 0.77). The value of B calcu- 
lated with many solutions remains zero to a lower temperature (T ~ 0.74) before it too 
becomes non zero, signalling either the onset of finite size effects or a breakdown of the 
approximation. In this temperature interval where deviations are clear yet the theory is 
reliable, the correlation length calculated from many solutions continues to grow while that 
of the maximal solutions rounds off. Certainly the rounded form coming from the maximal 
solutions is not a finite size effect as the curve for maximal solutions at 64"^ lies directly 
over this set. On the other hand the growing curve obtained from many solutions can 
become divergent as the size is increased since the number of solutions included will also 
increase. An exponent for this potential divergence can be determined by fitting the curve 
^ ~ (T — Tc)~'^- To do this, a combination of points from different data sets that include 
as many solutions as available, is used. The two lowest temperature points where B is 
non-zero are dropped, while at the other end of the range the fit is found to be insensitive 
to the inclusion or not of high temperature points. We find to be 1.25 ± 0.11 where the 
errors are statistical from jacknife. The value of Tc found in this way is 0.64 ± 0.01. 

The fit to the correlator, specifically the coefficient A, allows an estimate of the expo- 
nent f}. Scaling requires that the correlators G{x, ky = 0,kz = 0) behave as ^^~^G{x/^), 
leading in the case of the fitting function ( |5.4|) , to A ~ To avoid the uncertainty in 

Tc we plot log{A/^^) against log(^) in fig. ^ In contrast to the situation for the correlation 
length this logarithmic plot is curved and the high temperature points adversely affect the 
exponent. The line shown on the figure is obtained as the slope of the last points which 
are consistent with B = 0. This provides an estimate of = 0.89 with a statistical error 
of ±0.10. The systematic error associated with this measurement is harder to evaluate, 
but from the tendency of the lowest temperature points that were omitted, it is surely in 
a direction that would increase fj. 
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An alternative method of evaluating ry would be directly through the susceptibility: 
x'^'^-^Y.(^^)(Sj) = l3V{mf = (3v(^Y.'^-<) ■ (5.6) 

i,j \ i / 

After the average over disorder we find that the susceptibility is given in terms of the 
magnetisation fiuctuations x*^^^ = PVioa^ yielding a curve similar to the one obtained from 
the correlator fit. This curve does not however allow a determination of the exponent 
because in the critical low temperature region, m^ is contaminated with parts coming 
from non-zero magnetisation causing the slope to actually turn positive for the points not 
consistent with B = 0. 

5.2. Connected Correlator, 

The analysis of the connected correlator requires more delicacy. There are two reasons 
for this; firstly the extra FDT term and secondly a strong sample dependence on the ran- 
dom field. By evaluating the quantity a;^/3(a;^)^ for the fitting parameters, one can see that 
in contrast to the disconnected case, the distributions are distinctly non-gaussian. This 
observation is reinforced by an inspection of the susceptibilities for individual magnetic 
fields where one often observes points far removed from the mean. 

The fuU connected correlator is given by C^^) -C^^) +C'^^^'^\ The FDT part C^^^^\ 
defined in terms of gfj = dmf/dhj according to the plane-plane version of ( |2.(j| ), cannot 
be determined directly from the magnetisations. In fact a separate iteration is required to 
solve the equation for g°j: 

g,j = - m1){Dg,j + 5,j). (5.7) 

It would be prohibitive to perform this iteration for each site j of the lattice. However for 
a small number of cases we have done the iteration for one site fixed and note a strong 
dependence on the chosen site. We find that gij is more singular in position space than 
the part C^^) - C^^) arising from the magnetisation, and that asymptotically gij always 
goes to zero. This is consistent with our theoretical prejudice that a fit of the FDT part 
alone should be done with the function arising from a propagator with a single pole in 
momentum space. At higher temperatures the FDT piece is dominant in the connected 
correlator, at lower temperatures it can plausibly be neglected. In that case, an analysis of 
(7(2) _ (^(1) alone leads to a correlation length, at the lowest few temperature points, that 
is within the errors of the same quantity, fig. ^, calculated for the disconnected correlator. 
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To determine the exponent r] it is more convenient and accurate to abandon the fuU 
correlator and to work with the susceptibihty. 

;^(2) _ =f]v w^m'^m") - {J2 w^m'^f j (5.8) 

\ a ot / 

ij a i a 

The advantage of working with the susceptibihty is that (7° = J2j qTj can be evaluated by 



a single iteration of the equation obtained by summing (|5.7|) : 



g, = P{1 - m\) {Dg, + 1) . (5.9) 

The term x^^^'^^ is not very sensitive to the number of solutions employed to calculate it. 
In fact, in the temperature region of interest the FDT term is small compared with the 
difference term and is well approximated by the contribution from the maximal solutions 
alone. On the other hand, the magnetisation difference term does show some dependence 
on the number of solutions; the dependence is still small, but leads to different estimates 
for the exponent. Both terms, x'-^-* — x*-^-* and X^^^'^\ and their sum x^^'^ are shown in 
fig. H. The crossover to pure behavior is very clear and using only the points below the 
crossover temperature we can attempt to determine the exponent rj. The analog of the 
method used in the disconnected case, plotting x*-*^-*/^^ against correlation length, fails 
because the curve turns and begins to rise at low temperature in the same way as for x^^^ 
based on m^ We therefore resort to a direct determination of the susceptibility exponent 
7 defined by x*-*^-* ~ (T — Tc)~'^ ■ Because of the small number of points and the large 
errors the points are not weighted according to the size of their errors in this fit. We 
obtain 7 = 1.66 ± 0.53 with chequerboard data. For the increased number of solutions 
corresponding to chequerboard followed by two sets of random pertubations we find that 7 
is shifted to 7 = 1.97±0.91, where the increased size of the error is due the reduced number 
of random fields in this sample. Using the scaling relation 7 = z^(2 — r]), we determine the 
exponent 77 to be 77 = 0.7 ± 0.5 and 77 = 0.4 ± 0.8 for these respective data sets. 
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6. Conclusion 



By the direct procedure of solving the mean field equations for the RFIM we have anal- 
ysed the natural mean field theory obtained by weighting the solutions with Boltzmann-like 
factors. Our main conclusion is that this theory can describe critical behavior at the fer- 
romagnetic transition. The critical divergences arising from the sum over many solutions, 
each of which has regular behavior. We have investigated both the number and properties 
of the solutions that contribute. 

The exponents we have determined are, v = 1.25 ± 0.11, f] = 0.89 ± 0.10, r\ = 0.4 ± 0.8 
where the errors are statistical. The most accurate determinations to date are from the 



work of Rieger and Young who use Monte Carlo and finite size scaling . Our estimate of 
V is consistent with their value of = 1.4 ±0.2. Although our result for 77 lies close to their 
value of 77 = 0.6 ± 0.1, our errors are so large that this hardly provides a stringent test of 
our method. Given the systematic errors involved in our determination of 77, which would 
tend to increase it, there is no contradiction with their value of ry = 1.04 ±0.08. As is clear 
from the decay of the disconnected correlator at criticality, f\ must be greater than 1 for 
consistency of this analysis. The fact that it is so close to 1 has led to suspicions of a first 
order transition |jI^,|jT^. However we find it impossible to fit the specific heat, as measured 
by the derivative of the energy, to a divergence, and suspect a cusp form instead. There is 
thus the puzzle of zero latent heat already noted in [|14|. In order to see how this effect in 
particular, and the analysis in general, depends on the choice of probability distribution 
for the random fields, we have been studying the gaussian distribution. Preliminary results 
suggest that conclusions are unaffected. 

Only including those solutions that arise from chequerboard seeds gave good results 
for the disconnected correlator. Further accuracy would come from an increase in the 
number of random field samples rather than an increase in the number of solutions. In the 
connected case we observed stronger corrections coming from including solutions beyond 
chequerboard, however even there it was mainly the sample dependence that prevented 
greater accuracy. 
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Figure Captions 

Fig. 1. Total number of solutions found. Squares for chequerboard seeds alone, trian- 
gles for chequerboard followed by one set of random pertubations added to che- 
querboard, open triangles for chequerboard followed by two such sets of random 
pertubations. 

Fig. 2. W'^ = X]q;(^")^> the inverse of the effective number of solutions. Circles are 
for maximal solutions alone, squares for chequerboard, and triangles for chequer- 
board followed by (2 sets of) additional pertubations. 

Fig. 3. Local differences between maximal solutions at temperatures just below Tjisbi 
isosurfaces of difference m+i — rri-i = 0.1. 

Fig. 4. Number (top line) and weight (bottom line) of solutions shown as a histogram 
against magnetisation for temperatures decreasing towards the right. 

Fig. 5. Constant part of the disconnected correlator B, the squared magnetisation. Cir- 
cles indicate results for maximal solutions alone, squares are for chequerboard 
solutions. 

Fig. 6. Correlation Length ^, from the disconnected correlator. Circles indicate results 
for maximal solutions alone, squares are for chequerboard solutions. The line 
shows the best fit, ignoring the lowest two points which are not consistent with 
S = 0. 

Fig. 7. The parameter A/^^, from the disconnected correlator. Each point corresponds 
to the data set with most solutions available at that temperature. The last point 
is not included in the fit as S 0. 

Fig. 8. Connected Susceptibility, (x^^-* — X*^"*^^) P^rt shown with squares (chequerboard 
data), FDT piece x^^"^^^ shown with circles for the maximal solutions only (the 
errors are small), the sum x^^^ shown with open triangles. The line shows the 
best fit using only low temperature points. 
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